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Abstract: After a brief introduction into basic aspects of the formulation of lat- 
tice regularized QCD at finite temperature and density we discuss our current 
understanding of the QCD phase diagram at finite temperature. We present results 
from lattice calculations that emphasize the deconfining as well as chiral symmetry 
restoring features of the QCD transition, and discuss the thermodynamics of the 
high temperature phase. 



1 Introduction 



Almost immediately after the ground-breaking demonstration that the numer- 
ical analysis of lattice regularized quantum field theories ||l| can also provide quan- 
titative results on fundamental non-perturbative properties of QCD H it has been 
realized that this approach will also allow to study the QCD phase transition |^|^ 
and the equation of state of the quark-gluon plasma jsj. During the last 20 years 
we have learned a lot from lattice calculations about the phase structure of QCD 
at finite temperature. In fact, we do understand quite well the thermodynamics 
in the heavy quark mass limit of QCD, the pure SU{3) gauge theory, and even 
have calculated the critical temperature and the equation of state in this limit with 
an accuracy of a few percent. However, it is only now that we start to reach a 
level of accuracy in numerical calculations of QCD thermodynamics that allows to 
seriously consider quantitative studies of QCD with a realistic light quark mass 
spectrum. An important ingredient in the preparation of such calculations is the 
development of new regularization schemes in the fermion sector of the QCD La- 
grangian, which allow to reduce discretization errors and also improves the flavour 
symmetry of the lattice actions. The currently performed investigations of QCD 
thermodynamics provide first results with such improved actions and prepare the 
ground for calculations with a realistic light quark mass spectrum. 

The interest in analyzing the properties of QCD under extreme conditions is 
twofold. On the one hand it is the goal to reach a quantitative description of the 
behaviour of matter at high temperature and density. This does provide important 
input for a quantitative description of experimental signatures for the occurrence 
of a phase transition in heavy ion collisions and should also help to understand 
better the phase transitions that occurred during the early times of the evolution 
of the universe. Eventually it also may allow to answer the question whether a 
quark-gluon plasma can exist in the interior of dense neutron stars or did exist 
in early stages of supernova explosions. For this reason one would like to reach a 
quantitative understanding of the QCD equation of state, determine critical param- 
eters such as the critical temperature and the critical energy density and predict 
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the modification of basic liadron properties (masses, decay widtfis) witli temper- 
ature. On tlie other hand the analysis of a complicated quantum field theory like 
QCD at non-zero temperature can also help to improve our understanding of its 
non-perturbative properties at zero temperature. The introduction of external con- 
trol parameters (temperature, chemical potential) allows to observe the response 
of different observables to this and may provide a better understanding of their 
interdependence As such one would, for instance, like to clarify the role of con- 
finement and chiral symmetry breaking for the QCD phase transition. In which 
respect is the QCD phase transition deconfining and/or chiral symmetry restoring? 
In how far can the critical behaviour be described by intuitive pictures based on 
percolation, bag or resonance gas models which have been developed for the QCD 
transition? We will discuss these qualitative aspects of the QCD thermodynamics 
and also present results on basic questions concerning the equation of state and the 
critical temperature of the transition which ask for quantitative answers. 

In the next section we give a short introduction into the lattice formulation of 
QCD thermodynamics. In Section 3 we discuss the basic structure of the QCD phase 
diagram at finite temperature as it is known from lattice calculations. Section 4 is 
devoted to a discussion of basic thermodynamic observables which characterize the 
QCD transition to the plasma phase and we will identify general properties which 
show the deconfining and chiral symmetry restoring features of this transition. In 
Section 5 we comment on different length scales characterizing the QCD plasma and 
try to establish the temperature regime where lattice calculations may make contact 
with perturbative approaches. A description of recent results on the QCD equation 
of state and the critical temperature of the QCD transition which emphasizes the 
quark mass and flavour dependence of these quantities is given in Sections 6 and 
7, respectively. A brief discussion of the problems arising in lattice formulations of 
QCD at non-zero baryon number density or chemical potential is given in Section 8. 
Finally we give our conclusions in Section 9 and describe a specific set of improved 
gauge and fermion lattice actions in an Appendix. 

2 The Lattice Formulation of QCD Thermodynamics 
2.1 The basic steps from continuum to lattice ... 

Starting point for the discussion of the equilibrium thermodynamics of QCD on the 
lattice is the QCD partition function, which explicitly depends on the volume {V), 
the temperature (T) and the quark number chemical potential (/i). It is represented 
in terms of a Euclidean path integral over gauge (A^) and fermion {if>, ip) fields. 



where Av and ip, ip obey periodic and anti-periodic boundary conditions in Eu- 
clidean time, respectively. The Euclidean action Se = So + Sf contains a purely 
gluonic contribution (Sg) expressed in terms of the field strength tensor, F^,^ = 
d^Ai, — dvA^ — ig[A^, A^], and a fermion part (Sf), which couples the gauge and 
fermions field through the standard minimal substitution, 




(1) 



Se {V, T, fi) = Sg{V,T) + Sf{V, T, /i) 



(2) 
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l/T 

Sg(V,T)= I dxo I d^x^TrF^.F^. (3) 



V 
l/T 



/ = 1 



Sf{V, T, /i) = J dxo y > ^ V/ (7m [9m - iff^M] + 'W/ - M7o) V'/- (4) 

V 



Here are the different quark masses for the Uf different quark flavours and g 
denotes the QCD coupling constant. 

The path integral appearing in Eq. |l| is regularized by introducing a four di- 
mensional space-time lattice of size A'^^ x A'',- with a lattice spacing a. Volume and 
temperature are then related to the number of points in space and time directions, 
respectively, 

V = {N^ af , T"^ = iV, a . (5) 

While the discretization of the fermion sector, at least on the naive level, is straight- 
forwardly achieved by replacing derivatives by finite differences, the gauge sector 
is a bit more involved. Here we introduce link variables U^iix) which are associated 
with the link between two neighbouring sites of the lattice and describe the parallel 
transport of the field A from site x to x -\- jla, 

?7a:,M = Pexp( / dx^A^ix)] , (6) 

where P denotes the path ordering. The link variables Ufj,{x) are elements of the 
SU{3) colour group. A product of these link variables around an elementary pla- 
quette may be used to define an approximation to the gauge action, 



M^(i.i) = i_ iRe 



= ReTrf/„.^[/„+^,,[/t^,_^[/I,, 
= ^F^,F^,+Oia') . (7) 

A discretized version of the Euclidean gauge action, which reproduces the contin- 
uum version up to cut-off errors of order a^, thus is given by the Wilson action 



0, 



I3Sg = I3 J2 '^^'''^^^ =^ / 'i''^ + ' 

0<fJ,<v<3 

where we have introduced the gauge coupling /3 — 

As is well-known the naive discretization of the fermionic part of the action, 
which is obtained by introducing the simple finite difference scheme to discretize the 
derivative appearing in the fermion Lagrangian, i.e. d^il)f{x) = [ipn+fi — ^n-p.) f^a, 
does in the continuum limit not reproduce the particle content one started with. 
The massless lattice fermion propagator has poles not only at zero momentum but 
also at all other corners of the Brillouin zone and thus generates 16 rather than 
a single fermion species in the continuum limit. One thus faces a severe species 
doubling problem. The way out has been to either introduce an explicit chiral 
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symmetry breaking term, which is proportional to ad^Tpf{x) and thus vanishes 
in the continuum hmit (Wilson fermions Q), or to distribute the components of 
the fermion Dirac spinors over several lattice sites (staggered fermions) [Q]. The 
staggered fermion formulation does not eliminate the species doubling problem 
completely. One still gets four degenerate fermion species. However, it has the 
advantage that it preserves a continuous subgroup of the original global chiral 
symmetry. In the massless limit the chiral condensate thus still is an order parameter 
for the occurrence of a phase transition at finite temperature. 

Progress has been made in formulating lattice QCD also with chiral fermion 
actions which do avoid the species doubling and at the same time preserve the 
chiral symmetry of the QCD Lagrangian. This can, for instance, be achieved by 
introducing an extra fifth dimension pi. At present, however, very little has been 
done to study QCD thermodynamics on the lattice with these actions Much 
more is known on the QCD thermodynamics from calculations with Wilson and 
staggered fermions. We will in the following present results from both approaches. 
However, to be specific we will restrict ourselves here to a discussion of the staggered 
fermion formulation introduced by Kogut and Susskind Q. The fermion action can 
be written as 

^ = ^XnQ^mXm . (9) 

nm 

where the staggered fermion matrix Q'^^ is given by 
^ 3 

M=l 

+ ^(^n+6,m'^".o -5„,,„+6f^I,o) e"^ + ^nmm, . (10) 

Here we have introduced the chemical potential jl on the temporal links . As the 
fermion action is quadratic in the Grassmann valued quark fields x sjcid x we can 
integrate them out in the partition function and finally arrive at a representation 
of Z{V,T,fj,) on a 4-dimensional lattice of size A^^ x Nr, 



Z{N„,Nr,P,m.q,il) = 



y"[|d(7„,.(detQ^^(m„A))"^'/^e-''^'= . (11) 



We have made explicit the fact that the staggered fermion action does lead to four 
degenerate fermion flavours in the continuum limit, i.e. taking the continuum limit 
with the action given in Eqs. ^ and |lo| corresponds to = 4 in Eq. As the 
number of fermion species does appear only as an appropriate power of the fermion 
determinant, which is true also in the continuum limit, one also may choose n/ 7^ 4 
in Eq. |l^. This is the approach used to perform simulations with different number 
of flavours in the staggered fermion formulation. 

For /X = the fermion determinant appearing in Eq. |l^ is real and positive. 
Standard numerical techniques, which rely on a probability interpretation of the 
integrand in Eq. ^ thus can be applied. For fl ^ the determinant, however, 
becomes complex. Although the contribution of the imaginary part can easily be 
shown to be zero, as it should to give a real partition function, the real part is no 
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longer strictly positive. This sign problem so far still constitutes a major problem 
in the application of numerical techniques to studies of QCD at non-zero baryon 
number density or non-zero chemical potential. We therefore will restrict our dis- 
cussion of QCD thermodynamics mainly to the case p. = and will come back to 
the problems one faces for /i 7^ in Section 8. 

2.2 ... and back from the lattice to the continuum 

The lattice discretized QCD action discussed above reproduces the continuum ac- 
tion up to discretization errors of 0{a^). In order to perform the continuum limit 
at constant temperature, we will have to take the limit (a — > 0, Nt — > 00) with 
T = I/NtU fixed. In particular, for bulk thermodynamic observables like the pres- 
sure and energy density, which have dimension [T**] this limit is rather cumbersome. 
All lattice observables are dimensionless and are thus calculated in appropriate units 
of the lattice spacing a. As a consequence a calculation of, e.g., the pressure will 
provide pa* and thus yields a numerical result which decreases in magnitude like 
. Numerical calculation, however, are always based on the analysis of a finite 
set of suitably generated gauge field configurations and thus produce results which 
have a statistical error. It therefore rapidly becomes difficult to calculate bulk ther- 
modynamic quantities on lattices with large temporal extent Nt. For this reason 
it is of particular importance for finite temperature calculations to be able to use 
actions which have small discretization errors and thus allow to perform calcula- 
tions on lattices with moderate temporal extent. Such actions have been developed 
and successfully applied in thermodynamic calculations for the pure SU{3) gauge 
theory. In the fermion sector appropriate actions, which reduce cut-off effects in 
the high temperature ideal gas limit, so far have only been constructed for stag- 
gered fermions. As an example we describe a specific choice of improved gauge and 
staggered fermion actions in more detail in an Appendix. 

As mentioned above we have to perform the continuum limit in order to elim- 
inate lattice discretization errors and to arrive finally at quantitative predictions 
for the QCD thermodynamics. Eventually we thus have to analyze our observables 
on different size lattices and extrapolate our results to — > 00 at fixed tempera- 
ture. Unless we perform calculations at a well defined temperature, e.g. the critical 
temperature, we will have to determine the temperature scale from an additional 
(zero-temperature) calculation of an observable for which we know its physical value 
(in MeV) . This requires a calculation at the same value of the cut-off (same values 
of the bare couplings). Of course, we know such a quantity only for the physical 
case realized in nature, i. e. QCD with two light up and down quark flavours and a 
heavier strange quark. Nonetheless, we have good reason to believe that certain ob- 
servables are quite insensitive to changes in the quark masses, e.g. quenched hadron 
masses^ (mn) or the string tension (a) are believed to be suitable observables to 
set a physical scale even in the limit of infinite quark masses (pure 5*17(3) gauge 

^ A physical observable O is calculated on the lattice as dimensionless quantity, 
which we denote here by O. Quite often, however, we will also adopt the cus- 
tomary lattice notation, which explicitly specifies the cut-off dependence in the 
continuum limit, e.g. rfiH = iriHa or a = aa^. 
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theory). We thus may use calculations of these quantities to define a temperature 
scale, 

Tj^fa = \|^^^Nr or T /ruH = l/rfiHN-r . (12) 

In the pure SU (3) gauge theory as well as in the massless hmit the lattice spacing 
is controlled through /3, the only bare coupling appearing in the Euclidean action. 
Asymptotically a and (3 are then related through the leading order renormalization 
group equation, 

uAl ~ (66o//3)"*'^^''°e-''/'^*° , (13) 
where the two universal coefficients are given by. 




and is a scale parameter which unambiguously can be related to the scale 
parameter in other regularization schemes, e.g. to A-jjg. The continuum limit thus 
is reached with increasing /3. 

In the case of non-zero quark masses one has in addition to insure that the 
continuum limit is taken along a line of constant physics. This can be achieved by 
keeping a ratio of hadron masses, for instance the ratio of pseudo-scalar and vector 
meson masses, mps/mv , constant while varying the couplings {j3, rriq). In the limit 
/3 —f oo this requires a tuning of the bare quark masses such that — > 0. We also 
note that for small quark masses the vector meson mass mv approaches a constant, 
mv = TO,p+0('mq), while the pseudo-scalar is the Goldstone-particlc corresponding 
to the broken chiral symmetry of QCD (pion). Its mass is proportional to the 
square root of mq. In the following we will quite often quote results as a function 
of mps/mv which is just another way for quoting results for different values of the 
quark mass. 



3 The QCD Phase Diagram at Finite Temperature 

At vanishing baryon number density (or zero chemical potential) the properties of 
the QCD phase transition depend on the number of quark flavours and their masses. 
While it is a detailed quantitative question at which temperature the transition 
to the high temperature plasma phase occurs, we do expect that the nature of 
the transition, e.g. its order and details of the critical behaviour, are controlled by 
global symmetries of the QCD Lagrangian. Such symmetries only exist in the limits 
of either infinite or vanishing quark masses. For any non-zero, finite value of quark 
masses the global symmetries are explicitly broken. In fact, in the case of QCD the 
explicit symmetry breaking induced by the finite quark masses is very much similar 
to that induced by an external ferromagnetic field in spin models. We thus expect 
that a continuous phase transition, which may exist in the zero or infinite quark 
mass limit, will turn into a non-singular crossover behaviour for any finite value of 
the quark mass. First order transitions, on the other hand, may persist for some 
time before they end in a continuous transition. Whether a true phase transition 
exists in QCD with the physically realized spectrum of quark masses or whether 
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in this case the transition is just a (rapid) crossover, again becomes a quantitative 
question which we have to answer through direct numerical calculations. 

Our current understanding ol the qualitative aspects of QCD phase diagram is 
based on universality arguments for the symmetry breaking patterns in the heavy 
PH ] as well as the light quark mass regime ||l^,^. In the limit of infinitely heavy 
quarks, the pure 517(3) gauge theory, the large distance behaviour of the heavy 
quark free energy, Fqq, provides a unique distinction between confinement below 
Tc and deconfinement for T > Tc. On a lattice of size A^^ x Nr the heavy quark 
free energy^ can be calculated from the expectation value of the Polyakov loop 
correlation function 



exp 



T 



{TtL^TtLI) , rT = \x-y\Nr (15) 



where and Ly represent static quark and anti-quark sources located at the 
spatial points x and y, respectively. 



{xo,x) . (16) 



For large separations (r — *■ oo) the correlation function approaches \{L)\^, where 
(L) = N^^(Y2 Tria:) denotes the Polyakov loop expectation value, which therefore 
characterizes the behaviour of the heavy quark free energy at large distances and 
is an order parameter for deconfinement in the SU{3) gauge theory, 

/T\ (~ <^ confined phase, T < Tc 

^ ' I > ^ deconfined phase, T > Tc ' ^ ' 

The effective theory for the order parameter is a 3-dimensional spin model with 
global Z{3) symmetry. Universality arguments then suggest that the phase transi- 
tion is first order in the infinite quark mass limit |l^ . 

In the limit of vanishing quark masses the classical QCD Lagrangian is in- 
variant under chiral symmetry transformations; for Uf massless quark flavours the 
symmetry is 

Ua{1) X SULiuf) X SUninf). 

However, only the SUijif) flavour part of this symmetry is spontaneously broken in 
the vacuum, which gives rise to (n^ — 1) massless Goldstone particles, the pions. The 
axial Ua{^) only is a symmetry of the classical Lagrangian. It is explicitly broken 
due to quantum corrections in the QCD partition function, the axial anomaly, and 
therefore gets replaced by a discrete Z{nf) symmetry at low temperature. The basic 
observable which reflects the chiral properties of QCD is the chiral condensate. 



In the T ^ limit this is just the heavy quark potential; at non-zero temperature 
Fqq does, however, also include a contribution resulting from the overall change 
of entropy that arises from the presence of external quark and anti-quark sources. 
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In the hmit of vanishing quark masses the chiral condensate stays non zero as 
long as chiral symmetry is spontaneously broken. The chiral condensate thus is an 
obvious order parameter in the chiral limit, 

, , r > <4> symmetry broken phase, T < Tc 

\XX) I ^ Q ^ symmetric phase, T > Tc ' ^ ' 

For light quarks the global chiral symmetry is expected to control the critical 
behaviour of the QCD phase transition. In particular, the order of the transition is 
expected to depend on the number of light or massless flavours. The basic aspects of 
the n/-dependence of the phase diagram have been derived by Pisarski and Wilczek 
from an effective, 3-dimensional Lagrangian for the order parameter^ 

Ccff = -^Tr{d^^^d^$) - im2Tr(<pt<?) + ygi (Tr(<?t^))' 

2 

+ ^52Tr((<?^<?)^) +c(det<^ + det<?i^) , (20) 



with <!> = i,j = l,...,nf. jCcff has the same global symmetry as the QCD 

Lagrangian. A renormalization group analysis of this Lagrangian suggests that the 
transition is first order for n/ > 3 and second order for Uf — 2. The latter, however, 
is expected to hold only if the axial 17a (1) symmetry breaking, related to the det $ 
terms in Eq. |2^, does not become too weak at Tc so that the occurrence of a 
fluctuation induced first order transition would also become possible. 

This basic pattern has indeed been observed in lattice calculations. So far no 
indication for a discontinuous transition has been observed for n/ = 2. The tran- 
sition is found to be first order for n/ > 3. Moreover, the transition temperature 
is decreasing with increasing Uf and there are indications that chiral symmetry is 
already restored in the vacuum above a critical number of flavours [ p^ . 

The anticipated phase diagram of 3-flavour QCD at vanishing baryon number 
density is shown in Fig. |l| An interesting aspect of the phase diagram is the occur- 
rence of a second order transition line in the light quark mass regime, the boundary 
of the region of first order phase transitions. On this line the transition is controlled 
by an effective 3-dimensional theory with global Z{2) symmetry [Q, which is not 
a symmetry of the QCD Lagrangian. As this boundary lies in the light quark mass 
regime it may well be that this second order transition, for which neither the chiral 
condensate nor the Polyakov loop will be the order parameter, is equally important 
for the critical or crossover behaviour of QCD with a realistic quark mass spectrum 
as the nearby critical point in the chiral limit. In particular, we note that the critical 
exponent a is positive for the 3-d, Z{2) symmetric models whereas it is negative for 
the 0(4) model. A nearby Z{2) symmetric critical point in the QCD phase diagram 
will thus induce larger density fluctuations than would be expected in the vicinity 
of the chiral critical point. It therefore will be important to determine in detail the 
location of the physical point in the QCD phase diagram. 



^ It should be noted that this ansatz assumes that chiral symmetry is broken at 
low temperatures. Instanton model calculations suggest that the vacuum, in fact, 
is chirally symmetric already for n/ > 5 M]. 
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Td ~ 270 MeV 



m?fj* ~ 300 MeV 



Fig. 1. The QCD phase diagram of 3-flavour QCD with degenerate (u,d)-quark 
masses and a strange quark mass nis- 



4 Deconfinement versus Chiral Symmetry Restoration 

As outlined in the previous section the two properties of QCD, which explain the 
basic features of the observed spectrum of hadrons, are also of central importance 
for the structure of the QCD phase diagram at finite temperature - confinement and 
chiral symmetry breaking. While the former explains why we observe only colour- 
less states in the spectrum the latter describes the presence of light Goldstone 
particles, the pions. The confining property of QCD manifests itself in the long 
range behaviour of the heavy quark potential. At zero temperature the potential 
rises linearly at large distanceaj, Vqq{r) ~ ar, where a ~ (425 MeV)^ denotes the 
string tension, and forces the quarks and gluons to be confined to a hadronic bag. 
Chiral symmetry breaking leads to a non- vanishing quark anti-quark condensate, 
{qq) ~ (250 MeV)^ in the vacuum. Inside the hadron bag, however, the condensate 
vanishes. At high temperatures the individual hadronic bags are expected to merge 
to a single large bag, in which quarks and gluons can move freely. This bag pic- 
ture is closely related to percolation models for the QCD phase transition mn. It 
provides an intuitive argument for the occurrence of deconfinement and chiral sym- 
metry restoration. A priory it is, however, not evident that both non-perturbative 
properties have to get lost at the same temperature. It has been speculated that 
two distinct phase transitions leading to deconfinement at Td and chiral symmetry 
restoration at could occur in QCD General arguments about the scales 

* Here large distances actually refer to r ~ 1 fm. For larger distances the sponta- 
neous creation of quark anti-quark pairs from the vacuum leads to a breaking of 
the string, i.e. the potential tends to a constant value for r ^ oo (see Fig. 4). 
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involved^ suggest that Td < T^. Two distinct phase transitions indeed have been 
found in QCD related models like the SU{3) gauge theory with adjoint fermions 
|l^]. In QCD, however, there seems to be only one transition from the low tem- 
perature hadronic regime to the high temperature plasma phase. In fact, as can 
be seen from Fig. |l] there is a wide range of parameters (quark masses) for which 
the transition is not related to any singular behaviour in thermodynamic observ- 
ables; instead of a phase transition one observes just a rapid crossover behaviour. It 
thus is legitimate to ask which thermodynamic properties change when one moves 
from the low to the high temperature regime and to what extent these changes are 
related to deconfinement and/or chiral symmetry restoration. 

In the previous section we have introduced order parameters for deconfinement 
in the infinite quark mass limit, (L), and chiral symmetry restoration in the limit 
of vanishing quark masses, (V'V')- Related observables, which also signal a sudden 
change in the long distance behaviour of the heavy quark potential or the chiral 
condensate as function of temperature, are the corresponding susceptibilities, the 
Polyakov loop susceptibility (xl) and the chiral susceptibility (Xm), 



The behaviour of these observables is shown in Fig. g for the case of two flavour 
QCD with hght quarks. This clearly shows that the gauge coupling at which the 
different susceptibilities attain their maxima, or correspondingly the points of most 
rapid change in (L) and {'>ptp} coincide. Calculations of these observables for QCD 
with three degenerate quark flavours^ have been performed for a wide range of quark 
masses [^^. They confirm that the location of maxima in both susceptibilities are 
indeed strongly correlated. Within statistical accuracy they occur at the same tem- 
perature, although the height of these maxima is strongly quark mass dependent. 
This is shown in Fig. ^ For large {mps /mv>0.9) and small {mps /mv<0.3) quark 
masses the Polyakov loop and chiral susceptibility, respectively, show a strong vol- 
ume dependence, which is indicative for the presence of first order phase transitions 
in these corners of the phase diagram. Using zero temperature string tension cal- 
culations the pseudo-scalar meson masses have been estimated at which the first 
order transitions end in a second order transition. For better orientation in the 
phase diagram these estimates, which at present are not well established and cer- 
tainly are still subject to lattice artifacts (discretization errors, flavour symmetry 
breaking) , are shown in Fig. |^. As can be seen there is a broad range of quark (or 
meson) masses for which the QCD transition to the high temperature phase is a 
non-singular crossover. 

As expected the first order phase transition in the large quark mass regime is 
most clearly visible in the behaviour of the Polyakov loop susceptibility, i.e. the 
fiuctuation of the order parameter for the confinement-deconfinement transition in 
the pure gauge (m, — > oo) limit. Similarly, the transition in the chiral limit is most 
pronounced in the behaviour of the chiral condensate and its susceptibility. This 

^ The hadronic bag is larger than the constituent quark bag of a current quark 

surrounded by its gluon cloud. 
® This corresponds to calculations along the dotted, diagonal line in the phase 

diagram shown in Fig. H. 





d 



(21) 



Lattice QCD at High Temperature and Density 11 




m /T = 0.08 




m /T = 0.08 



Fig. 2. Deconfinement and chiral symmetry restoration in 2-flavour QCD: Shown 
is (L) (left), which is the order parameter for deconfinement in the pure gauge 
limit (niq — > oo), and (tpi^) (right), which is the order parameter for chiral sym- 
metry breaking in the chiral limit (m, ^ 0). Also shown are the corresponding 
susceptibilities as a function of the coupling /3 = 




Fig. 3. Quark mass dependence of the Polyakov loop and chiral susceptibilities 
versus mps/mv for 3-flavour QCD. Shown are results from calculations with the 
improved gauge and staggered fermion action discussed in the Appendix. 



emphasizes the chiral aspects of the QCD transition. One thus may wonder in what 
respect this transition in the light quark mass regime is a deconfining transition. 



4.1 Deconfinement 

When talking about deconfinement in QCD we have in mind that a large number of 
new degrees of freedom gets liberated at a (phase) transition temperature; quarks 
and gluons which at low temperature are confined in colourless hadrons and thus 
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Fig. 4. The left hand figure shows the heavy quark free energy in units of the 
square root of the string tension for the SU(3) gauge theory (open symbols) and 
three flavour QCD with light quarks (full symbols) . The right hand figure gives the 
hmiting value of the free energy normalized to the value at distance r ~ 0.23 fm 
p^ ] as a function of temperature for the case of three flavour QCD. The quark mass 
used in the n/ = 3 calculations corresponds to a ratio of pseudo-scalar and vector 
meson masses of mps/mv — 0.7. 



do not contribute to the thermodynamics, suddenly become liberated and start 
contributing to bulk thermodynamic observables like the energy density or pressure. 
In the heavy quark mass limit (rrig = oo) a more rigorous statement is based 
on the analysis of the long distance behaviour of the heavy quark free energy, 
which approaches a constant above Tc but diverges for T < Tc- For finite quark 
masses, however, the long distance behaviour of Fqq can no longer serve as an 
order parameter, the heavy quark free energy stays finite for all temperatures. 
This is shown in Fig. ^ where we compare results from calculations in the pure 
gauge theory jioj with results from a calculation in three flavour QCD [^. It is 
apparent from this figure that there is no drastic qualitative change in the structure 
of the heavy quark free energy as one crosses the transition temperature. Although 
the most rapid change in AV = Fqq{oo) — Fqq(r = O.S/y^) occurs for T ~ Tc 
and thus resembles the behaviour of the Polyakov loop expectation value and its 
susceptibility shown in Fig. |^, the heavy quark free energy does not seem to be a 
good indicator for deconfinement in the presence of light quarks. It looses its role 
as a rigorous order parameter and also does not refiect changes in the number of 
partonic degrees of freedom contributing to the thermodynamics. 

On the other hand, we expect that bulk thermodynamic quantities like the 
pressure do refiect the relevant number of degrees of freedom contributing to the 
thermodynamics in the high temperature limit. Due to asymptotic freedom the 
QCD pressure will approach the ideal gas value at infinite temperature. In this 
limit the number of degrees of freedom (quarks+gluons) is much larger than the 
three light pions which dominate the thermodynamics at low temperature. 
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Fig. 5. The pressure in QCD with different number of degrees of freedom as a func- 
tion of temperature. The curve labeled (2+l)-flavour corresponds to a calculation 
with two light and a four times heavier strange quark mass |2l|]. 



This change of active degrees of freedom is clearly visible in calculations of 
e.g. the pressure in the pure gauge sector and for QCD with different numbers of 
flavours. As can bee seen in Fig. B the pressure strongly reacts to changes in the 
number of degrees of freedom. It is this drastic change in the behaviour of the 
pressure or the energy density which indicates that the QCD (phase) transition to 
the plasma phase indeed is deconfining. However, it also is worthwhile to note that 
the transition does, in fact, take place at rather small values of the pressure (and 
energy density) . Only for temperatures T>2Tc does the pressure come close to the 
ideal gas limit so that one can, with some justification, identify the corresponding 
light degrees of freedom. This is the case for QCD with light quarks as well as in 
the quenched limit. At least for temperatures up to a few times Tc the dynamical 
degrees of freedom are certainly not just weakly interacting partons. 



4.2 Chiral symmetry restoration 

As chiral symmetry restoration does not lead to a significant change of light degrees 
of freedom, it also is not expected to have an appreciable effect on bulk thermo- 
dynamic observables - apart from controlling details of the transition very close to 
Tc- In particular, we expect that in the case of a continuous transition for Uf — 2, 
the chiral order parameter and its derivative, the chiral susceptibility, show critical 
behaviour which is characteristic for 0(4) spin models in three dimensions [^2|. The 
expected critical behaviour follows from standard scaling arguments derived from 
the singular part of the free energy density, 

fs(t,h) = -'^\nZs=b-''fs{b''H,b'"^h) , (23) 
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where t — {T — Tc)/Tc ~ {f3 — f3c) is the reduced temperature, h = m/T ~ rnqN^ the 
scaled quark mass and b is an arbitrary scale factor. For the chiral order parameter, 
and the chiral susceptibihty, Xm, one finds from Eq. 

{i)^) = h^^^F{z) (24) 

(25) 



F{z) - ^F\z) 



with scaling functions F and F' that only depend on a specific combination of the 
reduced temperature and scaled quark mass, z = th~^^'^^ . The critical exponents 
P and 5 are given in terms of yt and as /3 = (1 — yh)/yt and 5 = yh/{i ~ yh)- As 
the t-dependence enters in Xm{t, h) only through z one also deduces that the line of 
pseudo-critical couplings defined through the location of the maximum of Xmi^, h) 
at fixed h is described by a universal scaling function, 

t^h) = z, h'^^' . (26) 

Although there is ample evidence that the phase transition in 2-flavour QCD is con- 
tinuous in the chiral limit, the evidence for the expected 0(4) scaling is, at present, 
ambiguous. The behaviour of the pseudo-critical couplings is, in general, consistent 
with the expected scaling behaviour. The information on the magnetic equation of 
state, Eq. p4l however, seems to depend on the fermion discretization scheme used 
to analyze the critical behaviour. While calculations with Wilson fermions yield 
almost perfect agreement with the universal form of the 0(4) magnetic equation of 
state 1^^, significant deviations have been found in the case of staggered fermions 
lli^ j . The failure of the scaling analysis in the case of staggered fermions is surprising 
as this staggered fermion action has a global 0(2) symmetry even for finite values 
of the lattice cut-off and as the 0(4) and 0(2) magnetic equations of state are quite 
similar [Q. this suggests that finite size effects still play an irnportant role, which 
is supported by a recent finite size scaling analysis [ESl . In Fig. M we show the finite 
size scaling behaviour of the chiral condensate, which has been reanalyzed in ps]] . 
It is consistent with 0(4) (or 0(2)) scaling behaviour. This aspect, however, clearly 
needs further studies. 

The changes of the chiral condensate below Tc and chiral symmetry restora- 
tion at Tc will have a strong infiuence on the light hadron spectrum. At Tc the 
pseudo-scalar mesons (pions) will no longer be Goldstone particles, they turn into 
massive modes (quasi-particle excitations?) above Tc- Long distance correlations of 
the chiral condensate decay exponentially with a characteristic length scale propor- 
tional to the inverse scalar meson mass. A diverging chiral susceptibility at Tc thus 
indicates that the scalar meson mass vanishes at Tc- The mass splitting between 
parity partners thus will decrease when the symmetry breaking reduces and finally 
will become degenerate at Tc- 

As indicated above the modifications of the hadron spectrum are refiected by 
the temperature dependence of appropriately chosen susceptibilities, which are the 
space-time integral over hadronic correlation functions in a given quantum number 
channel, 

l-l/T I- 

XH = j dr J d'r GH{r,r) , (27) 
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Fig. 6. Finite size scaling of the chiral condensate in 2-flavour QCD |^^. Shown 
are data from calculations with standard staggered fermions on lattices of size 
8^ X 4 (circles), 12^ x 4 (triangles) and 16^ x 4 (squares). The calculations have 
been performed with different values of the quark mass at fixed value of the scaling 
variable Zc corresponding to the pseudo-critical line (Eq. p^ ). 

where the hadronic correlation function GH{T,r) for mesons is given by, 

GH{T,r) ^ {x{0)rHX{Omr,r)rHx{r,r)) , (28) 

and Fh is an appropriate combination of 7-matrices that projects onto a chosen 
quantum number channel. In particular, we note that the chiral susceptibility, Xm, 
defined in Eq. ^ is the susceptibility of the scalar correlation function. These sus- 
ceptibilities define generalized masses, mjj^ = xh, which are shown in Fig. |^ They, 
indeed, show the expected behaviour; scalar (/o) and pseudo-scalar (tt) partners be- 
come degenerate at Tc whereas the vector meson (5) which is related to the scalar 
meson through a C/a(1) rotation only gradually approaches the other masses. The 
axial (7a (1) symmetry thus remains broken at Tc. 

5 Screening at High Temperature — Short versus Long 
Distance Physics 

Our picture of the thermodynamics in the high temperature phase of QCD largely 
is influenced by perturbative concepts - asymptotic freedom and the screening of 
electric and magnetic components of the gluon fields. Asymptotic freedom suggests 
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Fig. 7. Temperature dependence of generalized hadron masses extracted from 
hadronic susceptibilities. Shown are results from calculations in 2-flavour QCD 
performed on lattices of size 8'^ x 4 with staggered fermions of mass = 0.02. 



that the temperature dependent running coupling, g{T), becomes small at high 
temperatures and eventually vanishes in the limit T — > oo. This in turn will lead 
to a separation of the thermal length scale, 1/T, from the electric, l/g(T)T, and 
mag netic, l/g^{T)T, screening length scales. The experience gained from lattice 
calculations in the pure SU (3) gauge theory, however, suggests that this separation 
of scales, unfortunately, will set in only at asymptotically large temperatures. For 
all interesting temperatures reachable in heavy ion experiments or even covering 
the temperature interval between the strong and electroweak phase transitions that 
occurred in the early universe, the coupling g{T) is of 0(1) and, moreover, the Debye 
screening mass is significantly larger than the leading order perturbative value. 



rriD 



1 + y g{T) T 



(29) 



In fact, for rc<r<100 Tc one finds that mo is still about three times larger than 
this leading order value j26| . A consequence of this large value of the screening mass 
is that also short distance properties of the plasma are strongly influenced by non- 
perturbative screening effects; the Debye screening length, ro = l/mo, becomes 
compatible with the characteristic length scale tsb in a free gas where the main 
contribution to the Stefan-Boltzmann law originates from particle with momenta 
p ~ ST, i.e. rsB ~ l/3r, and is of the same order as the mean separation be- 
tween partons in a quark-gluon plasma. We thus must expect that non-perturbative 
screening effects also have an influence on bulk thermodynamics properties (pres- 
sure, energy density) above Tc and even up to quite large temperatures. In fact, the 
calculations of the pressure and energy density in the SU{3) gauge theory and in 
QCD with light quarks, which we are going to discuss in the next section, show that 
at temperatures a few times Tc deviations from the ideal gas limit are still too large 
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to be understood in terms of conventional high temperature perturbation theory, 
which converges badly at these temperatures just because of the large contribution 
arising from Debye screening [Q. 

The screening of static quark and anti-quark sources is commonly analyzed in 
terms of Polyakov-loop correlation functions, which define the heavy quark free 
energy introduced in Eq. |l^. The leading perturbative contribution to Fqq{r,T) 
results from the exchange of two gluons, 

V,,{r,T) _ F,q(T,T)-F^q ^ {TrL,Tr4) 



T T \ \{L)\ 

1 (9\T) 



16 V Stt rT 



+ . (30) 



Higher order contributions will lead to screening of this powerlike large distance 
behaviour, 1/rT —> exp {—rrinr) /rT , and also result in an exponentiation of the 
leading order contribution. We then may split Fqq in contributions arising from 
quark anti-quark pairs in singlet (Fi) and octet (Fs) configurations 

g--F,-<j{'-.T)/T ^ 1 ^-Fi{r,T)/T ^ 8 ^-Fit(r,T)/T ^ j-g^j 

In accordance with zero temperature perturbation theory the singlet free energy 
is attractive whereas the octet free energy is repulsive. Their relative strength is 
such that it leads to a cancellation of the leading O(q ^) contributions to the colour 
averaged heavy quark free energy Fqq. From Eq. pll it is, however, apparent, that 
the cancellation of singlet and octet contributions only occurs at large distances. At 
short distances the contribution from the attractive singlet channel will dominate 
the heavy quark free energy, 

^iillll) = EliL:!! + eonst. 
T T 

= h const. for rT << 1 . (32) 

Stt rT 

In order to eliminate the subleading power-like behaviour at large distances we show 
in Fig. ^ (rT)'^ Vqq{r,T)/T calculated for the 5(7(3) gauge theory. As can be seen 
the change from the Coulomb-like behaviour at short distances to the exponential 
screening at large distances can be well localized. For Tc < T<2Tc it occurs already 
for rT ~ 0.2 or r ~ 0.15 {Tc/T) fm and shifts slightly to smaller rT with increasing 
temperature. A consequence of this efficient screening at short distances is that 
even heavy quark bound states get destroyed close to Tc in the plasma phase (J/i/"- 
suppression [p8|). 

The perturbative analysis of the heavy quark free energy also suggests that for 
fixed rT the only temperature dependence of V{r,T)/T arises from the running of 
the coupling g(T). The rapid change of V{r,T)/T at fixed rT which is apparent 
in Fig. H thus also suggests that for temperatures T<STc the coupling g(T) varies 
much more rapidly than the asymptotically expected logarithmic running with T. 

We thus conclude that non-perturbative screening effects are important for the 
thermodynamics in the plasma phase also for short distance observables which are 
sensitive to the physics at distances r>l/5T. The strong temperature dependence 
observed for r<3rc, moreover, suggests that the system cannot be described at a 
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Fig. 8. The heavy quark free energy at various temperatures in the deconfined 
phase of the SU (3) gauge theory. Calculations have been performed on lattices of 
size 32^ X 8 (filled symbols) and 32^ x 16 (open symbols) 



weakly coupled, asymptotically free plasma at these temperatures. These general 
features will carry over to the temperature dependence of the QCD equation of 
state which we are going to discuss in the next section. 



6 The QCD Equation of State 

The most fundamental quantity in equilibrium thermodynamics is, of course, the 
partition function itself, or the free energy density, 

f = -LinZ{T,V) . (33) 

All basic bulk thermodynamic observables can be derived from the free energy 
density. In the thermodynamic limit we obtain directly the pressure, p = —f and 
subsequently also other quantities like the energy (e) and entropy (s) densities or 
the velocity of sound (cs), 

e - 3p ^ _d_ f _p_\ _s_ ^ e+p „2 ^ dp 



In the limit of infinite temperature asymptotic freedom suggest that these observ- 
ables approach the ideal gas limit for a gas of free quarks and gluons, e — 3p = —3/ 
with p/T* given by Eq. M. Deviations from this ideal gas values have been studied 
in high temperature perturbation theory. However, it was well-known that this ex- 
pansion is no longer calculable perturbatively at 0{g^) |^9|. By now all calculable 
orders up to 0{g^ In g) have been calculated [Q. Unfortunately it turned out that 
the information gained from this expansion is rather limited. The expansion shows 
bad convergence behaviour and suggests that it is of use only at temperatures sev- 
eral orders of magnitude larger than the QCD transition temperature. In analytic 
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approaches one thus has to go beyond perturbation theory which currently is be- 
ing attempted by either using hard thermal loop resummation techniques [ ^ , p2[ 
or perturbative dimensional reduction combined with numerical simulations of the 
resulting effective 3-dimensional theory p^]. 

In order to make use of the basic thermodynamic relations, Eqs. ^ and 
in numerical calculations on the lattice we have to go through an additional in- 
termediate step. The free energy density itself is not directly accessible in Monte 
Carlo calculations; e.g. only expectation values can be calculated easily. One thus 
proceeds by calculating differences of the free energy density at two different tem- 
peratures. These are obtained by taking a suitable derivative of In Z followed by an 
integration, e.g. 



j_ 



This ansatz readily translates to the lattice. Taking derivatives with respect to the 
gauge coupling, /3 = rather than the temperature as was done in Eq. we 

obtain expectation values of the Euclidean action which can be integrated again to 
give the free energy density. 



J_ 

y4 



= iV.^ / d/3'({S)-{^)T=o) . (36) 

/3o "'/So 



Here 



is calculated on a lattice of size A''^ x A''t and (...)t=o denotes expectation values cal- 
culated on zero temperature lattices, which usually are approximated by symmetric 
lattices with Nt = N^. The lower integration limit is chosen at low temperatures 
so that //T* is small and may be ignored^. 

A little bit more involved is the calculation of the energy density as we have 
to take derivatives with respect to the temperature, T — l/N-ra. On lattices with 
fixed temporal extent Nr we rewrite this in terms of a derivative with respect to 
the lattice spacing a which in turn is controlled through the bare couplings of the 
QCD Lagrangian, a = a{f3,mq). We thus find for the case of n/ degenerate quark 
fiavours of mass ma 



rp4 T 



( d/9(a) 
y dlna 



(5) - {S)t=o 



7 



In the gluonic sector the relevant degrees of freedom at low temperature are 
glueballs. Even the lightest ones calculated on the lattice have large masses, 
ma — 1.5 GeV. The free energy density thus is exponentially suppressed already 
close to Tc. In QCD with light quarks the dominant contribution to the free 
energy density comes from pions. As long as we are dealing with massive quarks 
also this contribution gets suppressed exponentially. However, in the massless 
limit clearly some care has to be taken with the normalization of the free energy 
density. 
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An evaluation of the energy density thus e.g. requires the knowledge of two 13- 
functions. These may be determined by calculating two physical observables in 
lattice units for given values of /3 and m,; for instance, the string tension, cra^ and 
a ratio of hadron masses, mps/mv = m-nlmp. These quantities will have to be 
calculated at zero temperature which then also allows to determine a temperature 
scale in physical units as given in Eq.Jl3. This forms the basis for a calculation of 
the pressure as shown already in Fig. pT 

The numerical calculation of thermodynamic quantities is done on finite lattices 
with spatial extent Na and temporal extent Nt- In order to perform calculations 
close to the thermodynamic limit we want to use a large spatial extent of the 
lattice. In general it has been found that lattices with N^^'^N-r provide a good 
approximation to the infinite volume limit. In addition, we want to get close to the 
continuum limit in order to eliminate discretization errors. Taking the continuum 
limit at fixed temperature requires to perform the limit N-r — > oo. In order to per- 
form this limit in a controlled way we have to analyze in how far lattice calculations 
of bulk thermodynamic observables are influenced by the introduction of a finite 
lattice cut-off, i.e. we have to understand the systematic cut-off effects introduced 
through the non-zero lattice spacing. These cut-off effects are largest in the high 
(infinite) temperature limit which can be analyzed analytically in weak coupling 
lattice perturbation theory. We thus will discuss this limiting case first. 



6.1 High temperature limit of the QCD equation of state 

In the high temperature limit bulk thermodynamic observables are expected to 
approach their free gas values (Stefan-Boltzmann constants). In this limit cut-off 
effects in the pressure and in turn also in the energy density (esB = 3 psb) be- 
come most significant. Momenta of the order of the temperature, i.e. short distance 
properties, dominate the ideal gas behaviour. 

As discussed in Section 2 the most straightforward lattice representation of 
the QCD partition function in terms of the standard Wilson gauge and fermion 
actions as well as the staggered fermion action leads to a systematic 0{a?) cut-off 
dependence of physical observables. At finite temperature the temperature itself 
sets the scale for these 0{a^) effects, which thus give rise to 0{{aT)'^ = l/N^) 
deviations of e.g. the pressure from the continuum Stefan-Boltzmann value, 

+ ^ + 0(iVr*) . (39) 

One can eliminate these leading order cut-off effects by using improved actions 
which greatly reduces the cut-off dependence in the ideal gas limit. In the Appendix 
we discuss a specific set of improved gauge and fermion actions. In the gauge sector 
one may in addition to the standard Wilson plaquette term ((1 x l)-action) also 
include planar 6-link terms in the action. This is done in the O(a^) tree level 
improved (1 x 2)-action, which eliminates the leading order cut-off dependence 
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Fig. 9. Cut-off dependence of tlie ideal gas pressure for the 517(3) gauge theory 
(left) and several staggered fermion actions (right). These actions are defined in the 
Appendix. Cut-off effects for the Wilson fermion action are compatible with those 
of the standard staggered fermion action. 



completely. On lattices with temporal extent Nt one finds for the deviation of the 
gluonic part of the pressure jsij , 

PoW + + , (ixl)-action 

PG,SB )l + 0{Nr'^) , (1 X 2)-action 

A similar reduction of cut-off effects can be achieved in the fermion sector through 
the use of improved actions. So far, however, improved fermion actions, which re- 
duce or eliminate the leading order cut-off effects have only been constructed in the 
staggered fermion formulation. The Naik action |^^, which in addition to the ordi- 
nary one-link term in the staggered action also includes straight three-link terms, 
completely eliminates the 0(7V~^) errors on the tree level (ideal gas limit). The 
p4-action discussed in the Appendix does not eliminate this correction completely. 
It, however, reduces its contribution drastically over those present in the standard 
staggered actiorj^. 



MNt) 

Pf,sb 



1-H.57(^) +8.47(j^) +0{N-^) , 1-link standard 

staggered action (41) 
1 + 0.007 {^Y + 1.07 + 0{N-'') , p4-action 

as can be seen in Fig. ^. Moreover, the p4-action has the advantage that it improves 
the rotational symmetry of the fermion propagator which in turn also reduces vio- 
lations of rotational symmetry in the heavy quark potential. 

Using the tree level improved gauge action in combination with the improved 
staggered fermion action in numerical simulations at finite temperature it is possi- 
ble to perform calculations with small systematic cut-off errors already on lattices 
with small temporal extent, e.g. A'^^ = 4 or 6. In actual calculations performed with 



We quote here only an approximation to the Ar-dependence obtained from a fit 
in the interval 10 < A^ < 16. 
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Fig. 10. Pressure of the SU(3) gauge theory calculated on lattices with diflterent 
temporal extent and extrapolated to the continuum limit. Shown are results from 
calculations with the standard Wilson (1 x l)-action js^ and several improved 
actions [^,^, which are defined in the Appendix. The broad band shows the 
approximately self-consistent HTL calculation of |4l[. 



various actions in the pure gauge sector one finds that for temperatures T<5Tc 
the cut-off dependence of thermodynamic shows the pattern predicted by the in- 
finite temperature perturbative calculation. The absolute magnitude of the cut-off 
effects, however, is smaller by about a factor of two. This, of course, is reassuring 
for the numerical calculations performed with light quarks, where such a detailed 
systematic study of the cut-off dependence at present does not exist. 



6.2 Thermodynamics of the SU(3) gauge theory 

Before entering a discussion of bulk thermodynamics in two and three flavour QCD 
it is worthwhile to discuss some results on the equation of state in the heavy quark 
mass limit of QCD - the SU(3) gauge theory. In this case the temperature de- 
pendence of the pressure and energy density has been studied in great detail, cal- 
culations with the standard action and various improved actions p7| , ^ , ^ 
have been performed, the cut-off dependence has explicitly been analyzed through 
calculations on lattices with varying temporal extent A'',- and results have been ex- 
trapolated to the continuum limit. In Fig. ^ we show some results for the pressure 
obtained from such detailed analyzes with different actions |^,^,^. This figure 
shows the basic features of the temperature dependence of bulk thermodynamic 
quantities in QCD, which also carry over to the case of QCD with light quarks. 
The pressure stays small for almost all temperatures below Tc] this is expected, as 
the only degrees of freedom in the low temperature phase are giueballs which are 
rather heavy and thus lead to an exponential suppression of pressure and energy 
density at low temperature. Above Tc the pressure rises rapidly and reaches about 
70% of the asymptotic ideal gas value at T = 2 Tc. For even larger temperatures 




Fig. 11. Energy density, entropy density and pressure of tlie SU(3) gauge theory 
calculated on lattices with different temporal extent and extrapolated to the con- 
tinuum limit. The dashed band indicates the size of the latent heat gap in energy 
and entropy density 



the approach to this limiting value proceeds rather slowly. In fact, even at T ~ 4 Tc 
deviations from the ideal gas value are larger than 10%. This is too much to be 
understood in terms of weakly interacting gluons as they are described by ordi- 
nary high temperature perturbation theory |^^. Even at these high temperatures 
non-perturbative effects have to be taken into account which may be described in 
terms of interactions among quasi-particles Q,^. In Fig. |l^ we show the result 
of a self-consistent HTL resummation which leads to good agreement with 
the lattice calculations for r>3Tc. Other approaches [ ^Jic] ] reach a similarly good 
agreement in the high temperature regime. 

Compared to the pressure the energy density rises much more rapidly in the 
vicinity of Tc- In fact, as the transition is first order in the 5(7(3) gauge theory 
the energy density is discontinuous at Tc with a latent heat of about l.bT^ 
In Fig. [nl we show results for the energy density, entropy density and the pressure 
obtained from calculations with the Wilson action which have been extrapolated 
to the continuum limit |36|. The delayed rise of the pressure compared to that of 
the energy density has consequences for the velocity of sound in the QCD plasma 
in the vicinity of Tc. It is substantially smaller than in the high temperature ideal 
gas limit. 



6.3 Flavour dependence of the QCD equation of state 

As shown in Eq. ^ the pressure in QCD with light quarks can be calculated along 
the same line as in the pure gauge sector. Unlike in the pure gauge case it, however. 
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Fig. 12. The pressure in two flavour QCD calculated with unimproved gauge and 
staggered fermion actions (open symbols) |^^, RG-improved gauge and clover im- 
proved Wilson action (full symbols) and the p4-action (improved gauge and 
improved staggered fermions, see Appendix) (full line) |2l]. The dashed band esti- 
mates the results in the continuum limit as described in the text. The horizontal 
lines to the right and left show the Stefan-Boltzmann values for an ideal pion gas 
and a free quark-gluon gas, respectively. 



will be difficult to perform calculations on lattices with large temporal extent. In 
fact, at present all calculations of the equation of state are restricted to lattices 
with Nr = 4 and 6 [ pl|j4^ , p^ . The use of an improved fermion action thus seems to 
be even more important in this case. Of course, an additional problem arises from 
insufficient chiral properties of staggered and Wilson fermion actions. This will 
mainly be of importance in the low temperature phase and in the vicinity of the 
transition temperature. The continuum extrapolation thus will be more involved in 
the case of QCD with light quarks than in the pure gauge theory and we will have 
to perform calculations closer to the continuum limit. Nonetheless, in particular for 
small number of flavours, we may expect that the flavour symmetry breaking only 
has a small effect on the overall magnitude of bulk thermodynamic observables. 
After all, for n/ = 2, the pressure of an ideal massless pion gas contributes less 
than 10% of that of an ideal quark-gluon gas in the high temperature limit. For our 
discussion of bulk thermodynamic observables the main source for lattice artifacts 
thus still seems to arise from the short distance cut-off effects, which we have to 
control. Additional confidence in the numerical results can be gained by comparing 
simulations performed with different fermion actions. 

The importance of an improved lattice action, which leads to small cut-off errors 
at least in the high temperature ideal gas limit is apparent from Fig. ^ where we 
compare the results of a calculation of the pressure in 2-flavour QCD performed with 
unimproved gauge and staggered fermion actions and the RG-improved gauge 
and clover improved Wilson action with results obtained with the p4-action dis- 
cussed in the Appendix. At temperatures above T ~ 2 Tc these actions qualitatively 
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reproduce the cut-off effects calculated analytically in the infinite temperature limit 
(see Section 6.1). In particular, it is evident that also the Clover improved Wilson 
action leads to an overshooting of the continuum ideal gas limit. This is expected 
as the Clover term in the Wilson action does eliminate 0{ag^) cut-off effects but 
does not improve the high temperature ideal gas limit, which is 0(g^'). The clover 
improved Wilson action thus leads to the same large 0{a^) cut-off effects as the 
unimproved Wilson action. The infiuence of cut-off effects in bulk thermodynamic 
observables thus is similar in calculations with light quarks and in the SU{3) gauge 
theory. This observation may also help to estimate the cut-off effects still present 
in current calculations with light quarks. In particular, we know from the analysis 
performed in the pure gauge sector that in the interesting temperature regime of a 
few times Tc the cut-off dependence seems to be about a factor two smaller than 
calculated analytically in the infinite temperature limit; we may expect that this 
carries over to the case of QCD with light quarks. This is the basis for the estimated 
continuum extrapolation of the rif — 2 results shown as a dashed band in Fig. |l2| 

In Fig. P we have already shown result for the pressure calculated in QCD with 
different number of flavours. This figure clearly shows that the transition region 
shifts to smaller temperatures as the number of degrees of freedom is increased. 
Such a conclusion, of course, requires the determination of a temperature scale that 
is common to all QCD-like theories which have a particle content different from that 
realized in nature. We have determined this temperature scale by assuming that the 
string tension is flavour and quark mass independent. This assumption is supported 
by the observation that already in the heavy quark mass limit the string tension 
calculated in units of quenched hadron masses, e.g. mpj^fa = 1.81 (4) 0], is in 
good agreement with values required in QCD phenomenology, ^/a ~ 425 MeV. 

At high temperature the magnitude of p/T* clearly reflects the change in the 
number of light degrees of freedom present in the ideal gas limit. When we rescale 
the pressure by the corresponding ideal gas values it becomes, however, apparent 
that the overall pattern of the temperature dependence of p/T* is quite similar in 
all cases. This is shown in Fig. In particular, when one takes into account that 
a proper continuum extrapolation in QCD with light quarks is still missing this 
agreement achieved with improved staggered fermions is quite remarkable. 

We also note that the pressure at low temperature is enhanced in QCD with light 
quarks compared to the pure gauge case. This is an indication for the contribution 
of hadronic states, which are significantly lighter than the heavy glueballs of the 
SU(2>) gauge theory. This behaviour is even more clearly visible in the behaviour 
of the energy density which is shown in Fig. where we show results obtained 
with improved staggeredQand Wilson fermions. We note that these calculations 
yield consistent estimates for the energy density at Tc, 

^^^(%±2)Tt . (42) 

This estimate for ec/T^, which also is consistent with results obtained from 
calculations with a standard staggered fermion action , is an order of magnitude 
larger than the critical value on the hadronic side of the transition in the pure gauge 
theory (see Fig. |l^). It is, however, interesting to note that when we convert this 

^ The figure for staggered fermions is based on data from Ref. Here a contri- 
bution to e/T'* which is proportional to the bare quark mass and vanishes in the 
chiral limit is not taken into account. 
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Fig. 13. The pressure in units of tiie ideal gas pressure for tlie 5(7(3) gauge tlieory 
and QCD witli various number of flavours. The latter calculations have been per- 
formed on lattices with temporal extent Nt = 4 using the p4-action defined in the 
Appendix. Results are not yet extrapolated to the continuum limit. 



result for tc in physical units, [MeV/fm^], this difference gets to a large extent 
compensated by the shift in Tc to smaller values. When going from the infinite 
quark mass limit to the light quark mass regime the QCD transition thus seems 
to take place at compatible values of the energy density, Ec — (0.5 — l)GeV/fm''. 
The largest uncertainty on this number at present arises from uncertainties on the 
value of Tc (see next section). However, also the magnitude of tc/T* still has to be 
determined more accurately. Here two competing effects will be relevant. On the 
one hand we expect Ec/Tc to increase with decreasing quark masses, i.e. closer to 
the chiral limit. On the other hand, it is likely that finite volume effects are similar 
to those in the pure gauge sector, which suggests that tc/r* will still decrease closer 
to the thermodynamic limit, i.e. for Na — » oo. 

In the 2-flavour calculations performed with improved Wilson fermions the 
pressure and energy density have been calculated for several values of the quark 
mass, which corresponds to different ratios of the pseudo-scalar to vector meson 
mass raps/mv. The results show no significant quark mass dependence up to 
mps/mv — 0.9. This meson mass ratio corresponds to pseudo-scalar meson masses 
of about 1.5 GeV, which is somewhat larger than the ^-meson mass. This suggests 
that the corresponding quark mass is compatible with that of the strange quark. 
The approximate quark mass independence of the equation of state observed in the 
high temperature phase thus is consistent with our expectation that quark mass 
effects should become significant only when the quark masses get larger than the 
temperature. 
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Fig. 14. The energy density in QCD. The upper (lower) figure shows results from 
a calculation with improved staggered (Wilson ^ij) fermions on lattices with 
temporal extent N-r = 4 (A^t =4, 6). The staggered fermion calculations have been 
performed for a pseudo-scalar to vector meson mass ratio of mps/mv = 0.7. 



7 The Critical Temperature of the QCD Transition 

As discussed in Section 3 the transition to the high temperature phase is continuous 
and non-singular for a large range of quark masses. Nonetheless, for all quark masses 
this transition proceeds rather rapidly in a small temperature interval. A definite 
transition point thus can be identified, for instance through the location of peaks in 
the susceptibilities of the Polyakov loop or the chiral condensate defined in Eq. 
For a given value of the quark mass one thus determines pseudo-critical couplings, 
Ppcijriq), on a lattice with temporal extent A'^^- An additional calculation of an 
experimentally or phenomenologically known observable at zero temperature, e.g. 
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a hadron mass or the string tension, is still needed to determine the transition 
temperature from Eq. |l^. In the pure gauge theory the transition temperature, 
again has been analyzed in great detail and the influence of cut-off effects has been 
examined through calculations on different size lattices and with different actions. 
From this one finds for the critical temperature of the first order phase transition 
in the pure SU (3) gauge theory. 



(43) 



SU(3) gauge theory : Tc/V^ = 0.637 ± 0.005 
Tc = (271 ± 2) MeV 



Already the early calculations for the transition temperature with light quarks 
p^ , p7| indicated that the inclusion of light quarks leads to a significant decrease 
of the transition temperature. However, these early calculations, which have been 
performed with standard Wilson and staggered fermion actions, also led to 
significant discrepancies in the results for Tc as well as the order of the transition. 
These differences strongly diminished in the newer calculations which are based 
on improved Wilson fermions (Clover action) , domain wall fermions [|o[ 

as well as improved staggered fermions (p4-action) ]19|] . A compilation of these 
newer results is shown in Fig. ^ for various values of the quark masses. In order 
to compare calculations performed with different actions the results are presented 
in terms of a physical observable, the meson mass ratio mps/mv- In Fig. we 
show Tc/mv obtained for 2-fiavour QCD while Fi g. ISl b gives a comparison of 
results obtained with improved staggered fermions llOffor 2 and 3-fiavour QCD. 
Also shown there is a result for the case of (2+l)-fiavour QCD, i.e. for two light 
and one heavier quark flavour degree of freedom. Unfortunately the quark masses 
in this latter case are still too large to be compared directly with the situation 
realized in nature. We note however, that the results obtained so far suggest that 
the transition temperature in (2+l)-fiavour QCD is close to that of 2-fiavour QCD. 
The 3-flavour theory, on the other hand, leads to consistently smaller values of the 
critical temperature, Tc{nf — 2) — Tc{nf =3) ~ 20 MeV. The extrapolation of the 
transition temperatures to the chiral limit gave 



2 - flavour QCD : Tc = 



3 - flavour QCD : Tc = 



(171 ± 4) MeV, clover-improved Wilson 

fermions [48] 
(173 ± 8) MeV, improved staggered 

fermions |l^ 
(154 ± 8) MeV, improved staggered 

fermions [19] 



Here nip has been used to set the scale for Tc. Although the agreement between 
results obtained with Wilson and staggered fermions is striking, one should bear in 
mind that all these results have been obtained on lattice with temporal extent A',- = 
4, i.e. at rather large lattice spacing, a ~ 0.3 fm. Moreover, there are uncertainties 
involved in the ansatz used to extrapolate to the chiral limit. We thus estimate 
that the systematic error on the value of Tc/nip still is of similar magnitude as the 
purely statistical error quoted above. 

We note from Fig. |l5| that Tc/mv drops with increasing ratio mps/mv, i.e. 
with increasing quark mass. This may not be too surprising as my, of course, 
does not take on the physical p-meson mass value as long as mps/mv did not 
reach is physical value (vertical line in Fig. |l5|a). In fact, we know that Tc/mv 



Lattice QCD at High Temperature and Density 



29 




Fig. 15. Transition temperatures in units of mv- The upper figure shows a col- 
lection of results obtained for 2-flavour QCD with various fermion actions while 
in the lower figure we compare results obtained in 2 and 3-fiavour QCD with the 
p4-action described in the Appendix. All results are from simulations on lattices 
with temporal extent A^,- — 4. The large dot drawn for mps/mv = indicates the 
result of chiral extrapolations based on calculations with improved Wilson as 
well as improved staggered |^| fermions. The vertical line in the upper figure shows 
the location of the physical limit, mps = = 140 MeV. 

will approach zero for mps/mv = 1 as Tc will stay finite and take on the value 
calculated in the pure 5(7(3) gauge theory whereas mv will diverge in the heavy 
quark mass limit. Fig. ^ thus does not yet allow to quantify how Tc depends on 
the quark mass. A simple percolation picture for the QCD transition would suggest 
that Tc{mq) or better Tc{mps) will increase with increasing m,; with increasing 
ruq also the hadron masses increase and it becomes more difficult to excite the low 
lying hadronic states. It thus becomes more difficult to create a sufficiently high 
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Fig. 16. The transition temperature in 2 (filled squares) and 3 (circles) flavour 
QCD versus mpsi using an improved staggered fermion action (p4-action). Also 
shown are results for 2-flavour QCD obtained with the standard staggered fermion 
action (open squares). The dashed band indicates the uncertainty on Tc/^/o in the 
quenched limit. The straight line is the fit given in Eq. 0. 



particle/energy density in the hadronic phase that can trigger a phase (percolation) 
transition. Such a picture also follows from chiral model calculations [ pl| . 

As argued previously we should express Tc in units of an observable, which itself 
is not dependent on m,; the string tension (or also a quenched hadron mass) seems 
to be suitable for this purpose. In fact, this is what tacitly has been assumed when 
one converts the critical temperature of the SU(3) gauge theory T^j \fa ~ 0.63 into 
physical units as has also been done in Eq. 

To quantify the quark mass dependence of the transition temperature one may 
express Tc in units of y^. This ratio is shown in Fig. |l^ as a function of mpsj \fo- As 
can be seen the transition temperature starts deviating from the quenched values for 
'mps ^ (6 — ~ 2.5 GeV. We also note that the dependence of Tc on mps /\f^ 

is almost linear in the entire mass interval. Such a behaviour might, in fact, be 
expected for light quarks in the vicinity of a 2"'' order chiral transition where the 
dependence of the pseudo-critical temperature on the mass of the Goldstone-particle 
follows from the scaling relation, Eq. Gq, 

Tc{m^) - Tc{Q) ~ ml'"' . (44) 

For 2-flavour QCD the critical indices are expected to belong to the universality 
class of 3-d, 0(4) symmetric spin models and one thus would indeed expect 1/ pS — 
0.55. However, this clearly cannot be the origin of the quasi linear behaviour which 
is observed for rather large hadron masses and seems to be independent of n/. 
Moreover, unlike in chiral models the dependence of Tc on mps turns out to 
be rather weak. The line shown in Fig. ^ is a fit to the 3-flavour data, which gave 

^^'^ +0.04(1) fe") . (45) 



Lattice QCD at High Temperature and Density 



31 



It thus seems that the transition temperature does not react strongly on changes 
of the hghtest hadron masses. This favours the interpretation that the contributions 
of heavy resonance masses are equally important for the occurrence of the transition. 
In fact, this also can explain why the transition still sets in at quite low temperatures 
even when all hadron masses, including the pseudo-scalars, attain masses of the 
order of 1 GeV or more. Such an interpretation also is consistent with the weak 
quark mass dependence of the critical energy density we found from the analysis of 
the QCD equation of state in the previous section. 

For the quark masses currently used in lattice calculations a resonance gas model 
combined with a percolation criterion thus provides an appropriate to describe the 
thermodynamics close to Tc. It remains to be seen whether the role of the light 
meson sector becomes more dominant when we get closer to the chiral limit. 



8 Finite Density QCD 

Finite density calculations in QCD are affected by the well known sign pro blem, i. e. 
the fermion determinant appearing in the QCD partition function, Eq. hll becomes 
complex for non-zero values of the chemical potential n and thus prohibits the use 
of conventional numerical algorithms. The most detailed studies of this problem 
have so far been performed using the Glasgow algorithm which is based on a 
fugacity expansion of the grand canonical partition function at non-zero /i, 

aV 

ZgcWT,T,V)= V z^ZBiT,V), (46) 



where z = exp(^/r) is the fugacity and Zb are the canonical partition functions 
for fixed quark number _B; q = 3, 6 for one species of staggered or Wilson fermions, 
respectively. After introducing a complex chemical potential in Zqc the canonical 
partition functions can be obtained via a Fourier transformation 



1 Z"^" 

Zb{T,V)^—I d<l>e"^'' ZGcii(l>,T,V) (47) 







]^df/„,. 



.UBe-"'^'' , (48) 



with 

1 f^" 

aB = — d.^e''^^ (detQ^(m„i9i))"^/^ . (49) 
2^ J o 

The use of this ansatz for the calculation of canonical partition functions as 
expansion coefficients for Zqc has been discussed in [^,^ . A new approach has 
been suggested recently, which combines simulations with imaginary chemical 
potential with an analytic continuation based on the Ferrenberg-Swendsen multi- 
histogram method ||57|| . 
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Fig. 17. Polyakov loop expectation value (left) calculated on x 2 and the heavy 
quark potential (right) calculated on 16"' x 4 lattices in quenched QCD at zero and 
non-zero baryon number, B/3. 



One thus may evaluate the canonical partition functions as expectation values with 
respect to a trial partition function that can be handled numerically, for instance 
the partition function of the pure 5*17(3) gauge theory, 

aV 

Zgc{ij-/T,T,V) = Zsu(3) ^ 2"^{a-B)su(3) ■ (50) 

However, this approach so far did not overcome the severe numerical difficulties. 
Like other approaches it suffers from the problem that expectation values have to 
be calculated with respect to another ensemble so that the importance sampling 
which is at the heart of every numerical approach samples the wrong region of phase 
space and thus may become quite inefficient. 

It thus may be helpful to approach the finite density problems from another 
perspective. A reformulation of the original ansatz may lead to a representation of 
the partition function which, in the ideal case, would require the averaging over 
configurations with strictly positive weights only, or at least would lead to a strong 
reduction of configurations with negative weights. 

An alternative formulation of finite density QCD is given in terms of canonical 
rather than grand canonical partition functions js^], i.e. rather than introducing 
a non-zero chemical potential through which the number density is controlled one 
introduces directly a non-zero baryon number (or quark number B) through Eq. ^ 
from which the baryon number density on lattices of size A^^ x Nr is obtained as 
ub/T^ = ^(Nt/Nct)'^ . Also this formulation is by no means easy to use in general, 
i.e. for QCD with light quarks. In particular, it also still suffers from a sign problem. 
It, however, leads to a quite natural and useful formulation of the quenched limit of 
QCD at non-zero density which may be a good starting point for generalizing 
this approach to finite values of the quark mass. In the following we briefly outline 
the basic ideas of this approach. 

8.1 Quenched limit of finite density QCD 

It had been noticed early that the straightforward replacement of the fermion de- 
terminant by a constant does not lead to a meaningful static limit of QCD [m . In 
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fact, this simple replacement corresponds to the static limit fermion flavours carry- 
ing baryon number B and ~B, respectively [^. This should not be too surprising. 
When one starts with QCD at a non-zero baryon number and takes the limit of 
infinitely heavy quarks something should be left over from the determinant that 
represents the objects that carry the baryon number. In the canonical formulation 
this becomes obvious. For —> oo one ends up with a partition function, which 
for baryon number B/3 still includes the sum over products of B Polvakov loops, 
i.e. the static quark propagators which carry the baryon number [^^. This limit 
also has some analogy in the grand canonical formulation where the coupled limit 
mq,ii —> oo with exp(/i)/2mg kept fixed has been performed |6o|j6llP] . 

As the baryon number is carried by the rather heavy nucleons in the confined 
phase of QCD we may expect that it is quit reasonable to approximate them by 
static objects. This may already provide valuable insight into the thermodynamics 
of QCD at non-zero baryon number density already from quenched QCD. 

In the canonical approach simulations at non-zero B can be performed on rel- 
atively large lattices and the use of baryon number densities up to a few times 
nuclear matter density is possible | p^ , p2[ . The simulations performed so far in the 
static limit show the basic features expected at non-zero density. As can be seen 
from the behaviour of the Polyakov loop expectation value shown in Fig. ^ the 
transition region gets shifted to smaller temperatures (smaller coupling /3). The 
broadening of the transition region suggests a smooth crossover behaviour at non- 
zero density. However, in a canonical simulation it also may indicate the presence 
of a region of coexisting phases and thus would signal the existence of a 1"* order 
phase transition. This deserves further analysis. 

Even more interesting is the behaviour of the heavy quark free energy in the 
low temperature phase. As shown in the right frame of Fig. ^the free energy does 
get screened at non-vanishing number density. The influence of static quark sources 
on the heavy quark free energy is similar in magnitude to the screening (string 
breaking) seen in QCD simulations at finite temperature in the low temperature 
hadronic phase (see Fig. ^ . At non-zero baryon number density we thus may expect 
similarly strong medium effects as at finite temperature. 



9 Conclusions 

We have given a brief introduction into the lattice formulation of QCD thermody- 
namics and presented a few of the basic results on the QCD equation of state, the 
critical parameters for the transition to the QCD plasma phase and properties of 
this new phase of matter. 

The thermodynamics of the heavy quark mass limit is quite well under control; 
we know the equation of state and the transition temperature with an accuracy 
of a few percent. We now also have reached a first quantitative understanding of 
QCD with light quarks, which at present still corresponds to a world in which the 

This is a well known limit in statistical physics. When deriving the non-relativistic 
gas limit from a relativistic gas of particles with mass m, the rest mass is splitted 
off from the chemical potential, = jj.nr + fh, in order to cancel the corresponding 
rest mass term in the particle energies. On the lattice rh = ln(2mq) for large bare 
quark masses. 
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pion would have a mass of about (300-500) MeV. This stiU is too heavy to become 
sensitive to details of the physics of chiral symmetry breaking. Nonetheless, lattice 
calculations performed with different lattice fermion formulations start to produce a 
consistent picture for the quark mass dependence of the equation of state as well as 
the influence of the number of light flavours on the phase transition and they yield 
compatible results for the transition temperature. In these calculations we learn to 
control the systematic errors inherent to lattice calculations performed with a finite 
lattice cut-off and start getting control over the effects resulting from the explicit 
breaking of continuum symmetries in the fermion sector. With these experience at 
hand we soon will be able to study the thermodynamics of QCD with a realistic 
spectrum of light up and down quarks and a heavier strange quark. 
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Appendix: Improved Gauge and Fermion Actions 

Improved gauge actions 

When formulating a discretized version of QCD one has a great deal of freedom in 
choosing a lattice action. Different formulations may differ by sub-leading powers 
of the lattice cut-off, which vanish in the continuum limit. This has, for instance, 
been used by Symanzik to systematically improve scalar field theories josj and has 
then been applied to lattice regularized SU{N) gauge theories In addition 

to the elementary plaquette term appearing in the standard Wilson formulation of 
lattice QCD larger loops can be added to the action in such a way that the leading 
0{a^g'^) deviations from the continuum formulation are eliminated and corrections 
only start in 0{a'^g^ ,a?g^). A simple class of improved actions is, for instance, 
obtained by addingplanar loops of size (fe, /) to the standard Wilson action (one- 
plaquette action) |34|. The simplest extension of the Wilson one-plaquette action 
thus is to include an additional contribution from a planar six-link Wilson loop. 





+ 




) 









(51) 



Combining this six-link contribution with the four-link plaquette term in a suitable 
way one can eliminate the leading 0{a^) corrections and arrives at a formulation 
that reproduces the continuum action up to ©(a*) corrections at least on the clas- 
sical level at 0{g^), 



I3Sg = P 



"-1,1 ' ' n.iiu 



(52) 



0<M<i^<3 
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with ci,i = 5/3 and ci,2 = —1/6. We caU this action the tree-level improved 
(1 X 2)-action. It may be further improved perturbatively by eliminating the leading 
lattice cut-off effects also at 0{g^), i.e. aj cf'j + S^^'^j, or by introducing 
non-perturbative modifications. A well-studied gluon action with non-pert urbative 
corrections is the RG-improved action introduced by Y. Iwasaki pa. This RG- 
action also has the structure of Eq. |52| but with coefficients c^^ = 3.648 and 
c^2 = —0.662. Of course, this action will still lead to 0{a^) corrections in the ideal 
gas limit. The A'^^-dependence of cut-off effects resulting from these actions is shown 
in Fig. ^. In Section 6 we also show some results from calculations with a tadpole 
improved actions This non-perturbative improvement amounts to c\^f = ci.i 
and a replacement of ci,2 by Ci^2 = l/6uo(/3) where, 

In the ideal gas limit this action still has the same cut-off dependence as the tree- 
level improved (1x2) action. 



Improved staggered fermion actions 

When discussing the improvement of fermion actions there are at least two aspects 
one has to take into account. On the one hand one faces problems with cut-off 
effects similar to the pure gauge sector; on the tree level the standard Wilson and 
Kogut-Susskind discretization schemes introduce 0{a^) which will influence the 
short distance properties of physical observables. On the other hand also the global 
symmetries of the continuum Lagrangian are explicitly broken at non-zero lattice 
spacing. This influences the long distance properties of these actions, e.g. the light 
particle sector (Goldstone modes) of the lattice regularized theory. Both aspects 
are of importance for thermodynamic calculations. The latter problem certainly is 
of importance in the vicinity of the QCD phase transition while the former will 
show up when analyzing the high temperature limit of the equation of state. 

In the case of staggered fermions both problems have been addressed and 
schemes have been developed that lead to a reduction of cut-off effects at short 
distances, i.e. high temperature, and also allow to reduce the explicit flavour sym- 
metry breaking of the staggered discretization scheme. 

A particular form of improved action used in recent calculations of the Bielefeld 
group is a staggered fermion action, which in addition to the standard one- link term 
includes a set of bended three-link terms. 



Spimg) = Cl Sl-li„k,fat{uj) + C3 S-i-link + WT-q / ^ XxXx 

■ nn 



y X y 



+ UJ 



E 
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+1^9 X! ^^^^ • (^"^^ 

X 

Here ri^{x) = (— l)^o+- +^f ~i denotes the staggered fermion phase factors. Fur- 
thermore, we have made exphcit the dependence of the fermion action on different 
quark flavours q, and the corresponding bare quark masses ruq, and give an intuitive 
graphical representation of the action. The tree level coefficients cf and appear- 
ing in Sf have been fixed by demanding rotational invariance of the free quark 
propagator at 0{p'^) ( "p4-action" ) [^. In addition the 1-link term of the fermion 
action has been modified by introducing "fat" links |6^] with a weight lu = 0.2. The 
use of fat links does lead to a reduction of the fiavour symmetry breaking close to 
Tc and at the same time does not modify the good features of the p4-action at high 
temperature, i.e. it does not modify the cut-off effects at tree level and has little 
influence on the cut-off dependence of bulk thermodynamic observables at 0{g^) 
in the high temperature phase [^. Further details on the definition of the action 
are given in |69| . 

We refer to this action with a fat 1-link term combined with the tree level im- 
proved gauge action as the p4-action. The A^r-dependence of cut-off effects resulting 
from this action is shown in Fig. H. 
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